Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

hgvs: A Python package for manipulating sequence variants using HGVS nomenclature: 2018 Update.

The Human Genome Variation Society (HGVS) nomenclature guidelines encourage the accurate and standard description of DNA, RNA, and protein sequence variants in public variant databases and the scientific literature. Inconsistent application of the HGVS guidelines can lead to misinterpretation of variants in clinical settings. Reliable software tools are essential to ensure consistent application of the HGVS guidelines when reporting and interpreting variants. We present the hgvs Python package, a comprehensive tool for manipulating sequence variants according to the HGVS nomenclature guidelines. Distinguishing features of the hgvs package include: (1) parsing, formatting, validating, and normalizing variants on genome, transcript, and protein sequences; (2) projecting variants between aligned sequences, including those with gapped alignments; (3) flexible installation using remote or local data (fully local installations eliminate network dependencies); (4) extensive automated tests; and (5) open source development by a community from eight organizations worldwide. This report summarizes recent and significant updates to the hgvs package since its original release in 2014, and presents results of extensive validation using clinical relevant variants from ClinVar and HGMD

Antibacterial membrane attack by a pore-forming intestinal C-type lectin.

Human body-surface epithelia coexist in close association with complex bacterial communities and are protected by a variety of antibacterial proteins. C-type lectins of the RegIII family are bactericidal proteins that limit direct contact between bacteria and the intestinal epithelium and thus promote tolerance to the intestinal microbiota. RegIII lectins recognize their bacterial targets by binding peptidoglycan carbohydrate, but the mechanism by which they kill bacteria is unknown. Here we elucidate the mechanistic basis for RegIII bactericidal activity. We show that human RegIIIα (also known as HIP/PAP) binds membrane phospholipids and kills bacteria by forming a hexameric membrane-permeabilizing oligomeric pore. We derive a three-dimensional model of the RegIIIα pore by docking the RegIIIα crystal structure into a cryo-electron microscopic map of the pore complex, and show that the model accords with experimentally determined properties of the pore. Lipopolysaccharide inhibits RegIIIα pore-forming activity, explaining why RegIIIα is bactericidal for Gram-positive but not Gram-negative bacteria. Our findings identify C-type lectins as mediators of membrane attack in the mucosal immune system, and provide detailed insight into an antibacterial mechanism that promotes mutualism with the resident microbiota

Structure of the proton-gated urea channel from the gastric pathogen Helicobacter pylori

Half the world's population is chronically infected with Helicobacter pylori(1), causing gastritis, ulcers and increased incidence of gastric adenocarcinoma(2). Its proton-gated inner-membrane urea channel, HpUreI, is essential for survival in the acidic environment of the stomach(3). The channel is closed at neutral pH and opens at acidic pH to allow rapid urea access to cytoplasmic urease(4). Urease produces NH(3) and CO(2) that neutralize entering protons and thus buffer the periplasm to pH ∼6.1 even in gastric juice at pH <2.0. Here we report the structure of HpUreI, revealing six protomers assembled in a hexameric ring surrounding a central bilayer plug of ordered lipids. Each protomer encloses a channel formed by a twisted bundle of six transmembrane helices. The bundle defines a novel fold comprising a two-helix hairpin motif repeated three times around the central axis of the channel, without the inverted repeat of mammalian urea transporters. Both the channel and the protomer interface contain residues conserved in the AmiS/UreI superfamily, suggesting preservation of channel architecture and oligomeric state in this superfamily. Predominantly aromatic or aliphatic side chains line the entire channel and define two consecutive constriction sites in the middle of the channel. Mutation of Trp153 in the cytoplasmic constriction site to Ala or Phe reduces the selectivity for urea compared to thiourea, suggesting that solute interaction with Trp153 contributes specificity. The novel hexameric channel structure described here provides a new paradigm for permeation of urea and other small amide solutes in prokaryotes and archaea